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Abstract 

We present new spherically symmetric solutions of an SU(2) Einstein- Yang-Mills model coupled to a 
doublet of scalar fields. Sequences of asymptotically flat, Yang-Mills-boson star-type configurations are 
constructed numerically by considering an appropriate time-dependent ansatz for the complex scalar field 
and a static, purely magnetic SU(2)-Yang-Mills potential. Both nodeless as well as solutions with nodes of 
the scalar field and gauge potential are considered. We find that these solutions share many features with 
the "pure" boson stars. 



1 Introduction 



Oh! 
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Although there is still no direct evidence for the existence of scalar fields, there are many theoretical reasons 
' that these fields might play an important role in the evolution and the structure of the Universe. It is possible 
f~| , that a fraction of the bosonic dark matter has collapsed to form stellar type objects - boson stars. Boson 
stars (BS) are localised, static configurations of gravitationally bound zero temperature scalar particles, the 
complex scalar field possessing a harmonic time dependence. The study of BS started with the work of Kaup 
PP and Ruffini and Bonazzalo 0, who found asymptotically flat, spherically symmetric equilibrium solutions 
of the Einstein-Klein-Gordon equations. These configurations are "macroscopic quantum states" and are only 
prevented from collapsing gravitationally by the Heisenberg uncertainty principle. 

The BS share many features with their fermionic counterparts, presenting however many interesting differ- 
ences. For example, BS also exhibit a critical mass and critical particle number. Later work considered the 
self-interacting case j3] EI El or a non-minimal coupling of the scalar field to gravity [HJ. Boson stars in the 
presence of a dilaton or an axidilaton have also been studied by various authors [Jj, as well as boson- fermion 
stars jHJ- 

All these models have demonstrated the same characteristic: new interactions tend to increase the critical 
values of mass and particle number, although the particular values are very model dependent. The stability 
against perturbations around the equilibrium state has been discussed also by a number of authors |5|-|12|. An 
extensive review of the boson star properties is given in . |14) and more recently in |15) . 

Jetzer and van der Bij extended the BS model to include the coupling with a U(l) gauge group ^J. Here, 
the scalar field regularizes the central singularity of the Reissner-Nordstrom solution, which otherwise would 
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necessarily be present in the pure Einstein-Maxwell theory. A natural generalization of these charged BS 
configurations is to consider a larger gauge group. In this case, there are regular configurations even in the 
absence of a scalar field, as proven by the Bartnik-McKinnon (BM) static, time-independent solution of the 
coupled SU(2) Einstein- Yang-Mills (EYM) equations [T2| (see ^Sl f° r a general discussion of the properties of 
this type of solutions). 

In this paper, we extend the analysis of ^B] to include a complex doublet of scalar fields coupled to an SU(2) 
non-abelian gauge field 1 . In the static case, and for a symmetry breaking scalar potential, this system admits 
sphaleron solutions, describing the top of the potential-energy barrier separating gauge-inequivalent classical 
vacua ^Hl- The BM configurations are recovered in the limit of vanishing scalar field |22j . 

Here we present both numerical and analytical arguments for the existence of a new type of solution of 
the coupled EYM-scalar field equations, which combines the basic properties of both BS and BM models. For 
any value of the scalar field at the origin, the BS solutions can be generalized to include a BM particle inside. 
Different from the sphaleron case, however, the magnitude of the scalar field is zero at infinity. 

The paper is structured as follows: in the next Section we present the general framework and analyse the 
field equations and boundary conditions. In Section 3 we present our numerical results. We conclude with 
Section 4, where our results are summarized. 



2 General framework and equations of motion 
2.1 Basic ansatz 

Our study of the EYM-scalar field system is based upon the action 

1 1 

(1) 



S = I d xy/-g n 



1 -K- (Z? At $) t ( J D^$) - V($) - -F" F afa/ 



16ttG v m ' v ' y ' 4 



where, following the standard model, we take the scalar field $ to be a complex doublet. G is the gravitational 
constant and g m is the determinant of the metric tensor. 
Here is the usual gauge-covariant derivative 

D^ = d,& + g{A%a a )*, (2) 

while 

F% = %K - d»A% + geabcA^Al (3) 

is the SU(2) field strength tensor and g is the gauge coupling constant. 

In this letter we will restrict our analyse to the case V($) = /i 2 $^$, where \x is the scalar field mass, without 
including a scalar self-interaction term. As found in 0], although the inclusion of a A|$| 4 term drastically changes 
the value of the maximum mass and the corresponding critical central density of the boson star solutions, the 
qualitative features are essentially similar to the non-self interaction case. 

The action is invariant under a global phase rotation $ — > 5>e — %a which implies the existence of a conserved 
current 

J" = %ST (pv* f )4 - (D»$)&) . (4) 
This gives an associated conserved charge, namely, the number of scalar particles: 

N = J d a x^J t . (5) 

Since we assume spherical symmetry, it is convenient to use the metric in Schwarzschild-like coordinates: 

ds 2 = — - + r 2 (d9 2 + sin 2 6d<p 2 ) - cr 2 (r)B(r)dt 2 , (6) 

am 



1 Asymptotically anti de-Sitter q-stars interacting with a nonabelian field have been studied recently in 1231 . However, the 
properties of these notopological solitons are rather different as compared with the solutions discussed in this paper. 
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where B(r) = 1 — 2m(r) / r. m(r) may be interpreted as the total mass-energy within the radius r; its asymptotic 
value gives the total ADM mass of the solutions. 

Similar to the pure BS case, we can define a configuration radius 

R = If J ^ rJt ■ (7) 

We consider only static gauge fields, and for simplicity we assume A® = 0. The most general spherically 
symmetric SU(2) ansatz is parametrized by 

_ 1_- W(r) A W{r) (x . . o(r). - 

gr gr gr 



while for the scalar field we choose 



e- iM y>(r)+iK(r)(x a a a )](^j . (9) 



It is well known that the Ansatz for the matter fields is plagued with a residual gauge symmetry. Along 
with |24|. we fix the gauge by imposing the axial gauge 

x l Ai = => a(r) = . (10) 

In this paper we consider the consistent truncation of the full model W(r) — K(r) — 0, which corresponds 
to imposing the symmetry of the fields under the parity operator. Note that the original sphaleron solution 
of Klinkhamer and Manton as well as the gravitating generalizations (221 have been constructed in this 
reduced ansatz. However, solutions with W(r) ^ 0, K(r) ^ (so-called "bi-sphalerons" ) have been constructed 
in j2Hl ED ■ The general spherically symmetric configurations have also been used in [21] . 

2.2 Field equations and boundary conditions 

Within this ansatz, the classical equations of motion can be derived from the following 2-dimensional action: 



J [**G g 2 \ 2r 2 J \ V 2 



(ii) 



2 r v ' Bo 2 

where the prime denotes the derivative with respect to r. 

The field equations reduce to the following system of four non-linear differential equations: 

/ 1 / „ (I^ 2 -l) 2 \ , 1 N9 9 9 , 4> 2 r 2 Co 2 

rri = 4ttG — [BW' 2 + ( o2 ' )+Br 2 0' 2 + -tf(W-l) 2 + fi 2 r 2 (f ) 2 + ^— r - 
\g \ 2r z J 1 Ba A 

nGa [W 12 2l>2 cj ) 2 r 2 Cj 2 \ 

1 + ^V 2 + ^W > (12) 



g 2 Y a 2 B 2 
r A 2 



(aBW)' = ^L(W 2 -l) + -g 2 a4, 2 (W-l) 

2,~,2 



(aBr 2 (f)')' = ^{W-lY+^ar 2 - ^ " 



oB 



The regularity of the solution at the origin, the finiteness of the ADM mass and the requirement that 
the metric (O approaches the Minkowski metric for r — > oo lead to definite boundary conditions for the four 
functions to, <t, W, <p. As far as the metric functions are concerned we have to impose 

to(0) = , cr(oo) = 1 . (13) 

The remaining functions have to obey 

W(0) = 1, 0'(O) = , W(oo) = (-l) p , c/>(oo) = , (14) 
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(where p is a positive integer), which are the usual boundary conditions for the BS and BM configurations, 
respectively. Thus, for finite energy solutions, the term <D 2 </> 2 in forces cf> to tend to zero at infinity, which 
excludes a "boson star-sphaleron" hybrid model with a scalar field potential. 

The field equations imply the following behaviour for r — > in terms of three parameters (b, (f>Q, ctq): 



0(r) = o + y(M a V+O(r 3 ) , W(r) = 1 - br 2 + 0(r 4 ) , (15) 
m(r) = ^-L_ (66 2 + ^( M 2 - 3o)V 2 )) r 3 + 0(r 4 ) , „(r) = + + ^f^) r 2 + 0(r 3 ) . 

A similar analysis for r — > oo reveals that the leading order asymptotic behaviour of the metric functions is 
determined by the gauge field only, while the scalar field asymptotic form is similar to the pure BS case: 

0(r) = 4> ir A e-V^-^ r + ... , with A = -1 + M ^ = ^ , (16) 

Vv 2 - 



W(r) = (-If - + ..., m(r)=M = £+••■, a = 1 - 



To perform numerical computations and order-of-magnitude estimations, it is useful to have a new set of 
dimensionless variables. This is obtained by using the following rescaling 



r -> rV^G/g , cj>{r) -> (/)/V8^G , m(r) — > m(r) VAttG/cj . (17) 

The system of differential equations (|12fl with the above boundary conditions then depends on only two 
coupling constant, namely: 

a = n = — — ,u = w = — — 18) 

9 9 Mpi g g M Pt 

where Mpi — \fG is the Planck mass scale. 

Note that the factor to could be absorbed into the definition of the metric function a and this would change 
the boundary conditions for a, with cr(oo) ^ 1. The system of differential equations 112|) would then depend 
only on a. 

We have however kept two coupling constants, a and lo. In our numerical computations, we have used a 
and <^(0) as free parameters and obtained lo as function of 0(0). 



3 Numerical results 

The equations of motion l|12|) have been solved for a range of values of the scalar field at the origin, (j>(0), and 
several values of a. As expected, these solutions have many features in common with the pure BS and BM 
solutions; they also present new features that we will point out in the discussion. 

For each choice of 4>q, the field equations have a solution with the right asymptotics for r — ► oo only when 
(er(0), W"(Q) ~ b) in the expansion at the origin (|15fl take on certain values. The value of the shooting parameter 
b decreases with increasing <f>Q. For the same value of 0o, different values of Co corresponds to different numbers 
of nodes of the scalar field. 

In the following, we denote by p the number of nodes of the gauge field functions and by n the number of 
nodes of the scalar field function. For the considered configurations we could not find non-trivial nodeless gauge 
field function. Thus p = implies W(r) = 1. 



3.1 Known configurations: Bartnik-McKinnon and "pure boson" star solutions 

If we set the scalar field to zero, the corresponding equation is trivially satisfied and we are left with the 
spherically symmetric Einstein- Yang-Mills system, whose solutions are the BM configurations JJj . They consist 
of a discrete family of smooth solutions uniquely characterized by the number of nodes p of the function W(r), 
with p > 1. For these solutions the metric function B(r) possesses a minimum at some finite value of r. 
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Figure 1. The mass-parameter M, the particle number N, the value of the metric function a(r) at the origin, er(0), and 
the parameter u> are represented as a function of the value of the scalar field at the origin <^(0) for "pure" boson stars 
(p = 0, n = 0) and p = 1, n = 0, 1 Yang-Mills-boson star configurations. 

The metric function a(r) increases monotonically from r = (with < cr(0) < 1) to r = oo with er(oo) = 1. 
Boson-star solutions discussed in this paper can be regarded as deformations (i.e. with <p(r) ^ 0) of each solution 
of the BM sequence. 

In the case when the Yang-Mills field is set to its vacuum value the parameter a can further be rescaled in 
the radial variable, and the equations redce to those of the "pure" boson-star. The scalar field can be gradually 
deformed by setting (f)(0) = 4>o > and the corresponding solution can be constructed numerically The 
configurations can be characterized namely by their mass (m(oo)) the values 0(0), er(0) and B m (i.e. the minimal 
value of the metric function B(r)) and by the value w. 

The behaviour of M and N as a function of 4>(0) is well known in this case; the mass and particle number 
rise with increasing (f>(0) to the maximum values M max = 0.633 and N max — 0.653 (see Figure 1). Then M and 
N decrease, oscillate a bit and approach asymptotically a constant value. 
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Figure 2. The particle number N is represented as a function of the effective radius R for the same configurations as 
in Figure f . 

When the value of the scalar field at the origin 0(0) becomes large, the value cr(0) decreases and tends 
asymptotically to zero. However, a new phenomenon appears in this limit. Namely the metric function B(r) 
develops a second minimum between r = and the global minimum. This will be discussed in more detail in 
the context of YM-boson stars. 

Figure 1 also illustrates the fact that BS solutions exist only on a definite interval of u> and that on some 
subinterval (0.76 > u > 0.85) several solutions with the same w can exist. Here we have examined only the 
case when the scalar field possesses no nodes (i.e. n = 0), but solutions exist for which the scalar field function 
possesses one or more zeros, i.e. n^O. 

3.2 Yang-Mills-Boson-star: p = 1, 2, 3 

We have constructed (again by imposing a fixed value for 0(0) — 0o > 0) families of solutions with n = 0, 1, 
2, respectively 3 nodes of the scalar field function and p = 1, 2, 3 nodes of the gauge field function for varying 
a. We first discuss the solution for a generic value of the mass parameter a, namely a = 1. The results turn 
out to be qualitatively the same for other values of this parameter. A discussion of the domain of existence of 
solution in the parameter a is presented at the end of this section. 

The data characterizing the n — and the n = 1 solutions is shown in Figures 1,2. For the case n = 2 the 
pattern remains qualitatively the same. These graphs present the same quantities as in the case of the "pure" 
boson stars. 

The main difference is that the p-th BM solution is approached in the limit 0(0) — > 0. Accordingly, the 
ADM- mass does not vanish in this limit, but is equal to that of the corresponding BM solution. The particle 
number behaves like (f>(0) 2 and as a consequence the star radius stays finite. Excluding this latter point, note 
that the behaviour of M and N as a function of 0(0) is similar to that of the "pure" boson stars. 

We see also that the mass and the parameter u) depend only weakly on 0(0) (this holds true also for B m 
which is not shown in the Figure), while u(0) decreases monotonically for increasing 0(0). As in the case 
without a non-abelian field, we see that the particle number and the mass have their extrema, in particular 
their maximum, at the same value of the central density (characterized by 0(0)), decreasing monotonically to 
a constant asymptotic value. The behaviour of the particle number (or mass) in dependence on the radius is 
similar to the "pure" BS case (see Figure 2). In both cases the particle number versus radius curve shows a 
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Figure 3. The profiles of typical Yang-Mills-boson stars with a = 1 are plotted for several values of p, n. 

counterclockwise inspiralling at the critical points corresponding to a maximum of N . Another point which is 
worth pointing out is that the various solutions exist only on definite intervals of the parameter to and that lu is 
not a monotonic function of 0(0). This parameter converges to a finite value (depending on n,p) for 0(0) — > 0. 
Note that the values of lu for a fixed value of 0(0) are smaller in the case of Yang-Mills boson stars as compared 
to the "pure" boson stars and for the Yang-Mills boson stars themselves smaller for n — as compared to n = 1. 

Another feature which the Yang-Mills boson star shares with the "pure" boson star is the fact that, when 
the scalar field's value at the origin becomes large enough, the metric function B{r) develops a second local 
minimum, situated inside the global minimum. This is illustrated in Figure 3, where the profiles of the functions 
B(r), u(r), W(r), 0(r) are shown for three different combinations of p and n. Note that for p — 1 and n = 2 
the mentioned second minimum, which we associate with a shell-like structure of the solutions, appears. 

We remark here that our numerical analysis was stopped due to numerical difficulties related to the smallness 
of 0(0) and we can thus not make any conclusive statement about whether the second minimum becomes the 
global minimum if 0(0) is small enough. 

We have also managed to construct solutions with two nodes in the gauge field, i.e. p = 2 and n nodes in the 
scalar field. The profiles of the functions N(r), cr(r), W(r) and 0(r) for the p = 2 and n — 3 solution are shown 
in Figure 3. The qualitative behaviour of N and M is the same as for p = 1. The most sensitive parameter to 
the number of nodes p is the frequency lu. Choosing 0(0) rather small (e.g. 0(0) = 0.1) and concentrating on 
the n = solution, we observe that lu ~ 0.3 for p — 1 and lu ~ 0.9 for p = 2. This has important consequence 
on the domain of existence of solutions in the a-0(O)-planc. 

Let us finally discuss this domain of existence. The precise determination of this domain is a huge task, 
which we don't aim at in this letter, however, we can make some statements about the qualitative behaviour. 
We first assume that 0(0) is fixed (say < 0(0) < 1.5) and a is varied. Our numerical analysis suggests that 
the parameter lu decreases considerably when a decreases. A consequence of this is that the quantity a 2 — lu 2 , 
which appears in the asymptotic analysis (see (|16|l ) becomes imaginary for a < a cr such that exponentially 
localized solutions cease to exist for such values. For p = 1, we find numerically a cr w 0.16 for 0(0) = 0.1 and 
a cr ~ 0.12 for 0(0) = 0.5. For p = 2 the critical value of a turns out to be larger, e.g. a cr ~ 0.85 for 0(0) = 0.1. 

We could hardly find any relevant critical phenomenon limiting the solution pattern when increasing a. For 
instance we checked that both the metric and the scalar field stay perfectly regular up to a ~ 10.0. 

For fixed a and increasing 0(0) we have found that the value of the metric function a at the origin, <r(0), 
becomes very small for large 0(0), while the minimum of the function B{r) stays finite. This suggests that 
the solution becomes singular at the origin for large enough 0(0). However, because the numerical calculations 
become very difficult for large 0(0), we refrain from making any further conclusions here. 
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4 Conclusions and further remarks 



In this paper we have investigated a new type of configuration combining the basic properties of two well 
known solutions: boson stars for a time-dependent complex scalar field coupled to gravity on the one hand and 
Bartnik-McKinnon particle-like solutions of the SU(2) Einstein- Yang-Mills system on the other hand. 

The properties of these Yang-Mills boson star solutions differ considerably from their abelian counterparts 
discussed in [T2|. Note also that these solutions represent the simplest examples of Yang-Mills boson star 
configurations. We expect that static, respectively rotating axially symmetric solutions as well as dilaton 
generalizations may exist as well. 

The Einstein- Yang-Mills system possesses also black hole solutions, so-called "coloured" black holes, which 
are counterexamples to the no-hair conjecture JHj. There are no black hole analogues to the "pure" BS con- 
figurations [2HJ though. However, the introduction of an SU(2) gauge field invalidates the arguments presented 
in |25| . The existence of black hole counterparts of the configurations discussed in this paper is thus an open 
question. 

An important physical question when discussing self-gravitating configurations is whether these solutions are 
stable. The stability of the Yang-Mills boson stars can be studied in general by considering a perturbation about 
the classical solutions and solving the linearized equations. This leads to a system of coupled Schrodinger (or 
Sturm-Liouvillc) equations, whose analysis is a complicated task. However several useful informations about the 
negative modes of the solutions can be obtained by arguments based of catastrophe theory and the inspection 
of bifurcations and/or cusps occuring in the pattern of the solutions |2lfl[T2] . 

It is known that the p-node BM solution is plagued with 2p unstable modes [57], in other words it is a 
sphaleron for p > 0. As a consequence it is very likely that, at least for small values of 0(0), the non-abelian 
boson star obtained by deformation of the p-th BM solution also posseses 2p unstable modes. (Note also that 
we could not find non-trivial p — solutions.) This is at least consistent with the stability of the boson star, 
i.e. for p = 0. 

According to the arguments of the catastrophe theory, we have to look for bifurcations or the occurence of 
cusps in the plot of e.g. the energy (or better the binding energy given by the difference M — aN) as a function 
of the particle number N. Here the main difference between "pure" boson stars and the Yang-Mills boson stars 
resides in the fact that for the latter the mass of the solution does not tend to zero in the limit N — > because 
it converges to the energy of the corresponding BM solution. As a consequence, the function M — aN is strictly 
positive for N << 1 suggesting that the binding energy is positive and contorting the statement that the p > 1 
Yang-Mills boson stars are indeed unstable. 

In the case p = 1 such a diagram indicates the occurence of a cusp for the maximal value of the particle 
number , i.e. at N w 0.35, corresponding to (f)(0) » 0.4 (the results here correspond to a = 1). No solution 
seems to exist for N > 0.35, however, the solutions existing for higher values of cf>(0), i.e. 0(0) > 0.4, form 
another branch with a cusp at N « 0.35. The binding energy of the second branch is higher than the one 
corresponding to the main branch. So it is likely that the solution constructed has two unstable modes for the 
first branch, i.e. for 0(0) < 0.4 and three unstable modes for 0(0) > 0.4. These results definitely need to be 
confirmed by a normal mode analysis of the linearized equations. 
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